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Abstract 

The morphology of Bucky onions is investigated by ab initio calcu- 
lations using Yang's new 0{N) method. It is found that for large 
single shell fullerenes with Ih symmetry, the spherical morphology 
has the lower energy than that of polyhedrons. The formation en- 
ergy per atom follows a simple scaling law. Including an estimate 
of inter-shell van de Waals interactions leads to the conclusion that 
spherical multiple shell Bucky onions are likely the most stable 
structure of large carbon clusters. These results are in good agree- 
ment with recent observation. 
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The exciting discovery of concentric spherical graphitic shells by Ugartea raises several 
intriguing questions. Most fundamentally, it suggests that multi-shell fuUerenes — so called 
Bucky onions — are the most stable structures of finite carbon clusters. This, and the 
remarkable spherical shape of these onions call for theoretical studies. For cluster size up to 

= 240 carbon atoms, Adams et a/.i have performed first principle calculations and showed 
that "ball shaped" fuUerenes have a lower energy than that of tube-shaped. However, based 
either on empirical potentials calculations, geometric consideration and elastic theory, it has 
been concluded by several groups that for large fuUerenes the morphology is polyhedrally 
facetedi^. This appears to contradict the experimental observation of spherical shells. 

In this letter, we report results of first principle calculations for carbon clusters up to 
1,000 atoms by the "divide- and-conquer" method of Yangi. These are the first ab initio 
calculations for such large molecules. The results unequivocally demonstrate that for single 
shell fuUerenes the spherical morphology has a lower energy than those of polyhedrons. 
Furthermore we show that the leading dependence of the total energy obtained can be 
understood in term of the Hiickel theory. Including the inter-shell van de Waals interactions 
leads to the conclusion that the multi-shell concentric spherical Bucky onions are likley the 
most stable structure of large carbon clusters, as observed in recent experiment Jil. 

Most methods of first principle calculations are limited to molecules of small or modest 
size. Recently one of us has developed an efficient 0{N) algorithm^, designed specifically 
for large molecules. The method has been implemented for general molecular computations 
and tested in small systems against the traditional Kohn-Sham density-functional theory. 
The tests performed were molecular bonding energetic, molecular internal rotation potential 
barriers, and total electronic density of stated. In our present calculations, a cluster is 
divided into subsystems with one atom each, and the local basis set for each subsystem used 
includes the atomic orbitals of the atom as well as those up to the third nearest neighbors. 
We use the local density approximation(LDA) for the exchange and correlation energy and 
numerical LDA solutions of a spherical carbon atom for the atomic orbitals. To reduce the 
computational effort the non-selfconsistent Harris functional! is employed as well as the Ih 



symmetry. Using this new method we are able to carry out ab initio calculations for carbon 
clusters up to 1, 000 atoms. Selective results are summarized in Table I. 

We restrict our investigation to Goldberg type I fullerene£l, in which the number of atoms 
in nth shell is N{n) = 60n^ and the average radius can be approximated by Rn ~ 2.46n, 
where b is the average C-C bond length. Thus, the inter-shell spacing is close to the inter- 
layer distance in the graphite, giving rise to the maximum inter-shell van de Waals attraction. 
For the smallest member Ceo, all sites are equivalent, the shape can be classified equally 
well as spherical or polyhedral. The next shell C240 has three independent atomic sites 
(Fig.l). The Ih symmetry reduces the independent variables to seven. A global search for 
the minimum is out of reach of our present computational capabilities. We have investigated 
four most likely morphologies. 1) Sphere (S), in which all three sites are assumed to be on 
a spherical surface. 2) Icosahedrally (I) faceted, here all three sites are assumes to be on a 
faceted surface. A generalization of this morphology to large fuUerenes is an Icosahedroni. 
3) Truncated Icosahedron (TI), this is an inflated version of Cgo, in which site 1 and 2 (Fig.l) 
are assumed to be on the pentagonal facet while 2 and 3 are on a hexagonal facet. 4) The 
low energy morphology obtained by Yoshida and Osawa (YS)i. This morphology is neither 
spherical nor simply faceted. However, as measured by the deviation from a sphere it is very 
close to (TI). 

For the morphologies (S), (I) and (TI) the total energy was minimized with respect to 
two independent bond lengths (Fig.l). For the (YS) shape we take the ratios of various 
bond lengths to be the same as that given by Ref.4. Table I lists the formation energy 
and corresponding geometric parameters obtained. As one can see, the formation energy 
for S-C240 is significantly smaller than that of I-C240, TI-C240 and YO-C240- Therefore, our 
results demonstrate that for C240 the spherical morphology has a lower energy than that of 
polyhedrons. This conclusion remains unchanged when the calculations were repeated using 
the atomic orbitals from up to fourth nearest neighbor atoms. 

Included in Table I are also results of similar investigations for C540. Again it is found 
that the spherical morphology has significant lower energy. Strong evidence in favor of 



the spherical shape is the fact that in all cases the formation energy is correlated with the 
sphericity as measured by the standard deviation (SD) from a perfect sphere (Table I). The 
more spherical is the morphology the lower is the energy. Based on these results we conclude 
that for large single-shell fullerenes the spherical shape is the preferred morphology. 

To understand the dependence of the formation energy several more clusters are 
examined. Plotted in Fig. 2 are the formation energy per atom Ec versus Ec decreases 

monotonically with increasing fullerene size N. A straight line fitting to the four largest 
fullerenes calculated leads to 

E,(iV) = -7.178(1.0-^) eV. (1) 

Thus, Ec{oc) = — 7.178eV is the extrapolated cohesive energy per atom for single graphitic 
layer. This number is close to that found by other calculations and by experimental. 

The simple dependence can be understood in terms of the Hiickel theory. In the 
nearest neighbor (n.n.) tight binding calculation there are two contributions to the energy 
difference between a fullerene and that of a graphitic layer. The first is due to the finite size 
of the fullerene. Exact calculation shows that the tight binding energy increases linearly 
with Et{N) = — 1.574t(l — ^^^^), where t is the n.n. hoping integralS. This relation 
should reflects the approximate dependence of a bonding energy. The second, and more 
important contribution comes from the fact that for vr orbitals the n.n. overlaps depends on 
the molecular sizei. For spherical morphology t^^ is proportional to the cos{(j)), where is the 
radial angle spanned by the nearest neighbor C-C bond, t^, = to cos = to cos(2 arcsin(2^)) 
where is the radius of the sphere and to is the limiting value for the graphitic plane. 
Taking Rj^ ^ 2Ah\J N /QQ and expanding in to the leading order one gets 

5.21 

t. ^ to(1.0 - ^) (2) 

Compared with Eq.l one sees the dominant contribution to the dependence of the cohesive 
energy is due to the change in orbital overlaps. 

We now discuss the energetics of multi-shell Bucky onions and their morphology. Here 
one needs to include the inter-shell interactions. In Goldberg type I fullerenes the spacing 
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between successive shells, ~ 3.4A, is almost the same as that of interlayer spacing in graphite. 
Therefore one expects that interactions are dominated by the van de Waals force. In solid 
Cqo, the intermolecular interactions have been successfully modeled using the van de Waals 
interactionJli. There it was found that 90% of cohesion energy can be calculated by simply 
summing over all the inter-fullerenes C-C interactions using the Lennard- Jones potential. 
For Bucky onions, this method should be an even better approximation as the inter-shell 
spacing are uniform. Therefore we use the same Lennard- Jones potential (e = 2.964 meV, 
a = 3.407A) to evaluate the inter-shell interactions. 

The fact that Rn+i — Rn ~ cr suggests that the inter-shell attraction is the largest when 
both shells have the spherical shape. Indeed numerical calculations show that the attrac- 
tive interactions in the sphere morphology is significantly larger than that of polyhedrons, 
because the non-uniform spacing between shells in the latter. For example, the attraction 
between C240 and S-C540 (-17.7 eV, see Table II) is significantly larger than that between 
I-C*24o aiid I-C540 (-12.3 eV). Thus, the inter-shell interactions also favor the spherical shape. 
Therefore we conclude that the shape of a multi-shell Bucky onion is spherical. This is in 
agreement with experimental observation of spherical shaped multi-shell concentric Bucky 
onionsB. In Fig. 3 a non-perspective view of a five shell Bucky onion is shown. 

From Table II the van de Waals interactions between the nth and {n+ l)th shell is found 
to be well described by Vn,n+i = — 3.0r;,(?T, -|- 1) Using the extrapolated energy Eq.l 

for the single shell fullerene, one can estimates that when n > 3 two shell Bucky onion 
with N{n) + N{n + 1) atoms is more stable than a single shell fullerene with the same total 
number of atoms. The first case where the number N{n) +N{n+ 1) also belongs to a type I 
Goldberg fullerene is when n = 3. Thus our calculations suggest that energetically C(54o,96o) 
is more stable than Cisoo^- This result is rather different from earlier calculation where the 
crossover was found to occur at ~ 600oi. Finally, it is clear that for large clusters the 
dense Bucky onions are most stable. This is because the van de Waals attraction is ~ for 
a / shell onion, while the cost is ~ / compared with a single shell fullerene. 
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In conclusion, we have shown that the shape of large fuUerenes are spherical regardless 
whether it is single or multiple shell. By including the van de Waals interactions between 
shells it is demonstrated that spherical dense Bucky onions are likely the most stable struc- 
ture for large carbon clusters. This is in good agreement with recent experimental observa- 
tions of these concentric spherical shells. The possibility of forming a lattice of Bucky onions 
should open up a new area of fuUerene research and many interesting physical properties 
can be expected. 
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true for type II in which N{n) = 20n^. 

^ In the case of spherical morphology, it is found that the rotational barrier of inner shell vs 
outer shell is very small. Thus, free rotation is expected at room temperature. In contrast, 
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^ The next occurrence for the identity n"^ + {n + 1)^ = is when n = 20, m = 29. In this 
case N — 50460, the van de Waals attraction 1^20,21 = 1261 eV is large and the two-shell 
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FIGURES 

FIG. 1. Asymmetric part of large fullerene structures, (a) C240, shown are three independent 
sites and two bond labeling (see text), (b) Similar drawing for C540. 

FIG. 2. The formation energy vs 1/A''. The straight line (Eq.l) is the least square fit into the 
four points = 180, 240, 540, 960. 

FIG. 3. A non-perspective view of the five shells dense Bucky onion C(6o,240,540,960,i500) • I^^r 
clarity, only portion of each shell is shown. 
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TABLES 



TABLE I. 


Formation energies per 


atom for various fullerenes with different morphologies: S - 


spherical, I 


Icosahcdron faceted, TI 


- truncated Icosahedron, YO 


the structure given 


in Ref.4. 


b - average bond length, R - Average 


radius, SD - standard deviation from a perfect sphere. 




Ec (eV) 


6(A) (61, 62) 


R{A) 


SD (A) 


S-C240 


-7.05 


1.43 (L43, 1.43) 


7.12 


0.00 


I-C'240 


-6.90 


1.45 (1.46, 1.42) 


7.27 


0.39 


TI-C240 


-6.92 


1.46 (1.47, 1.43) 


7.09 


0.18 


YO-C240 


-6.97 


1.45 (1.45, 1.40) 


7.03 


0.17 


S-C540 


-7.09 


1.41 (1.42, 1.41) 


10.5 


0.00 


I-C540 


-6.86 


1.43 (1.44, 1.40) 


10.4 


0.52 


TABLE IL The van de Waals interactions between the first five shells of a dense 


spherical 


Bucky onion. 










Vij (eV) 


C240 


C540 


C960 


C1500 


Ceo 


-5.57 


-0.81 


-0.20 


-0.07 


C240 




-17.7 


-2.34 


-0.55 


C540 






-35.9 


-4.28 


C960 








-59.9 
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